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Abstract 

The three moments of inertia associated with the wobbhng mode built on the superdeformed 
states in ^^^Lu are investigated by means of the cranked shell model plus random phase approxi- 
mation to the configuration with an aligned quasiparticle. The result indicates that it is crucial to 
take into account the direct contribution to the moments of inertia from the aligned quasiparticle 
so as to realize Jx > Jy in positive-gamma shapes. Quenching of the pairing gap cooperates with 
the alignment effect. The peculiarity of the recently observed ^^'^Lu data is discussed by calculating 
not only the electromagnetic properties but also the excitation spectra. 
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Rotation is one of the specific collective motions in finite many-body systems. Most of 
the nuclear rotational spectra can be understood as the outcome of one-dimensional (ID) 
rotations of axially symmetric nuclei. Two representative models — the moment of inertia of 
the irrotational fluid, JT'"'^'^, and that of the rigid rotor, JT"^, both specified by an appropriate 
axially-symmetric deformation parameter (3 — could not reproduce the experimental ones, 
^cxp. jiTT ^ jj-cxp ^ From a microscopic viewpoint, the moment of inertia can 

be calculated as the response of the many-body system to an externally forced rotation 
— the cranking model This reproduces J'^^^ well by taking into account the pairing 
correlation. Triaxial nuclei can rotate about their three principal axes and the corresponding 
three moments of inertia depend on their shapes in general. In spite of a lot of theoretical 
studies, their shape (in particular the triaxiality parameter 7) dependence has not been 
understood well because of the lack of decisive experimental data. Recently some evidences 
of three-dimensional (3D) rotations have been observed, such as the shears bands and the 
so-called chiral-twin bands [0. In addition to these fully 3D motions, from the general 
argument of symmetry breaking, there must be a low-lying collective mode associated with 
the symmetry reduction from a ID rotating axially symmetric mean field to a 3D rotating 
triaxial one. This is called the wobbling mode. Notice that the collective mode associated 
with the "phase transition" from an axially symmetric to a triaxial mean field in the non- 
rotating case is the well known gamma vibration. Therefore the wobbling mode can be 
said to be produced by an interplay of triaxiality and rotation. The wobbling mode is 
described as a small amplitude fluctuation of the rotational axis away from the principal 
axis with the largest moment of inertia. Bohr and Mottelson first discussed this mode [^. 
Mikhailov and Janssen and Marshalek ||^ described this mode in terms of the random 
phase approximation (RPA) in the rotating frame. In these works it was shown that at 
7 = this mode turns into the odd-spin members of the gamma-vibrational band while 
at 7 = 60° or —120° it becomes the precession mode built on top of the high-K isomeric 
states 0. Here we note that, according to the direction of the rotational axis relative to the 
three principal axes of the shape, 7 runs from —120° to 60°. 

Recently electromagnetic (EM) properties of the second triaxial superdeformed (TSD2) 
band in ^^^Lu were reported and concluded that the TSD2 is a wobbling band excited on 
the previously known yrast TSDl band, on the basis of comparisons to a particle- rotor 
model (PRM) calculation [0, ^. In conventional PRM calculations an irrotational moment 
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of inertia, 

J-r = ^:7osin^(7 + ^vrfc), (1) 

where k = 1-3 denote the x, y and z principal axes, is assumed. The magnitude J'o is treated 
as an adjustable parameter although it can be identified as J'o = 8-82/?^, where B2 is the 
inertia parameter in the Bohr Hamiltonian f^]. This reduces to J'"'^ in the first paragraph 
by substituting 7 = and k = 1, and satisfies such a required property that collective 
rotations about the symmetry axes are forbidden. Since Jy" is largest for < 7 < 60° and 
the main rotation occurs about the axis of the largest inertia, the PRM with J^" can not 
describe the positive-gamma rotation, that is, the rotation about the shortest axis (x axis). 
Then in Refs. [|^, |^ the so-called gamma-reversed moment of inertia |Ty], JT^*^^, defined by 
inverting the sign of 7 in Eq.(|I]), was adopted. Although this reproduced the measured EM 
properties well, this does not satisfy the required property mentioned above and its physical 
implications are not very clear. In this Letter, therefore, we study the moments of inertia 
associated with the wobbling motion excited on the positive-gamma states by means of the 
cranked shell model plus RPA. This framework does not divide the system into a valence 
particle and a rotor, and therefore can calculate the three moments of inertia of the whole 
system microscopically. We believe that this is the first step toward understanding the fully 
3D nuclear rotations. 

We have developed a computer code for the RPA to excitation modes built on configura- 
tions with arbitrary number of aligned quasiparticles (QPs). In this Letter we present the 
results for the 4 - 6QP configurations in Gd isotopes and the IQP one in ^^^Lu. In particu- 
lar, this is the first RPA calculation for the rotating odd-A configurations, to our knowledge. 
Note that this approach is different from the conventional particle-vibration coupling cal- 
culations where the RPA itself is performed for the even-even "core" configurations. Since 
the details of the formulation have already been given in Refs. JTT|, O, here we describe 



only the outline. The QP states were obtained by diagonalizing the cranked triaxial Nilsson 
plus BCS Hamiltonian at each rotational frequency cUrot with adjusting chemical potentials 
to give correct average particle numbers. The doubly-stretched F and 1 ■ s potentials were 
adopted, and their strengths were taken from Ref . . The RPA calculation was performed 
with adopting the pairing plus doubly-stretched Q ■ Q interaction. The existence of aligned 
QPs is taken into account by exchanging the definition of the QP creation and annihila- 
tion operators in an appropriate manner. Actual calculations were done in five major shells 



(^j\[iosc) _ 2 _ -J? gj^^ N^"^^^ = 2 - 6) by using the dispersion equation p, 

{hu) = {hu,,,) IreffU .^rnffW . ' (2) 



obtained by decoupling the Nambu-Goldstone mode analytically assuming 7 7^ 0. This 
equation is independent of the strengths of the interaction. Not only the collective wobbling 
mode {oj = u^oh) but also many non-collective modes are obtained from this equation. 
The effective inertia J'y^Pi^^) = J^p\uj) /Vty^z{'^), defined in the principal-axis frame (their 
concrete expressions were given in Ref. |jl2|]), depend on the eigenmode while the kinematical 
<Jx = {Jx)/^rot: where the expectation value is taken with respect to the whole system, is 
common to all modes. It should be noted that Eq.(^ coincides with the original expression 
for uj^oh if Jx and Jy'^^^ (1^) are replaced with constant moments of inertia. 

In the following, we present some numerical results. Here the parameters €2 (alternative 
to /?), 7, An and Ap were chosen so as to reproduce available experimental data, and kept 
constant as functions of cUrot- We have confirmed that qualitative features of the result are 
robust and the details of the parameter dependence will be given in a separate publica- 



tion [|I4|. It is non-trivial to obtain the wobbling solution in the RPA for positive-gamma 
nuclei and the QP alignment is indispensable for its appearance. In order to show this, we 
first discuss a theoretical calculation for a precession mode which might be built on top of 
the I'^ = 49/2"*" isomeric state in ^^''Gd, where the whole angular momentum is built up by 
the alignment of the five QPs, [(7r/;,ii/2)^(z//i9/2, /7/2)^]i8+ in ^"^^Gd plus [1^^13/2] 13/2+ , so that 
a 7 = 60° shape (axially symmetric about the x axis) is realized. This state is obtained 
by cranking with /itUrot = 0.3 MeV. We chose 62 = 0.19 and An = Ap = 0.6 MeV, and 



reproduced the observed static quadrupole moment and (^-factor |T^, |T^. In order to see 
the behavior of the three moments of inertia, we calculated the wobbling mode by changing 
the parameter 7 from 60°. The result is presented in Fig.|I|(a). Although at a first glance 
their 7 dependence resembles that of the rigid rotor, 

='-^^,(^1-^^ cos {^ + lnk)y (3) 

the physical contents of J'x changes with 7; the fraction of the collective contribution de- 
creases as 7 increases and reaches at 7 = 60°. Accordingly it can be conjectured that the 7 
dependence of the "rotor" contribution is approximately irrotational and the QP contribu- 
tion is superimposed on top of the former by aligning its angular momentum to the x axis. 



Our previous calculation |T2|, ^ for a negative-gamma nucleus, ^^^Os, also supports this and 
consequently it is thought that the wobbling mode can appear relatively easily in superfiuid 
negative-gamma nuclei. To see if this conjecture is meaningful, starting from ^^^Gd we add 
the ^13/2 quasineutrons sequentially. The result shows that increases as the number of 
aligned QPs increases. Since the increase of Jy^f^ is rather moderate, the increase of 
leads to that of the wobbling frequency uj^^h- Thus, the exchange from < Jy in 
to Jx > Jy in may be related qualitatively to the increase of Jx stemming from the 
alignment which is not accounted for in the PRM, considering the fact that the alignment 
of particle states leads to 7 > 0. 

At 7 ~ 30° where Jy reaches its maximum as in the irrotational model, we could not 
obtain a wobbling solution. In Fig.|l|(b), Tiu^ob and the wobbling angle, 

Owoh = tan — , (4) 

are graphed. This shows that w^oh becomes imaginary and 9^oh blows up in this region. 
Comparing Fig.|l](a) and (b), it may be inferred that the wobbling motion excited on a mean 
field rotating about the x axis becomes unstable at 7 ~ 30° due to J^ < Jy^^\ and that a 
tilted-axis rotation would be realized. Putting this unstable region in-between, the solution 
in the larger-7 side is like a precession of an axially symmetric body about the x axis, 
whereas that in the smaller-7 side is like a gamma vibration around an axially symmetric 
shape about the 2; axis. 

Now we turn to the TSD bands in ^^^^Lu. We chose 62 = 0.43, 7 = 20° and = Ap = 
0.3 MeV, and obtained transition quadrupole moments Qt = 10.9 - 11.3 eb for hu^ot = 0.20 
- 0.57 MeV in accordance with the data, Qt = 10.7 ± 0.7 eb [0. We have obtained for the 
first time (aside from the theoretical simulation above) the wobbling solution in the RPA for 
positive-gamma nuclei. Here it should be stressed that the inclusion of the five major shells 
and the alignment effect of the proton ^13/2 quasiparticle is essential for obtaining this result. 
In Fig.^(a) the measured excitation energy of the TSD2 band relative to that of the TSDl 
and the calculated hu^oh are shown. The most peculiar point in the experimental data is 
that cUwob decreases as a function of cUrot- If t^rot-independent moments of inertia such as 
the irrotational ones are adopted, Wwob increases linearly with c^rot, see the comment below 
Eq.(^. The wobbling frequency is sensitive to the difference among the three moments of 
inertia, and the ratios J^^^^ Jx and J^'^^y Jx actually determine a;„ob- For example, the 



gamma- reversed moments of inertia give Jy""^ j = 0.43 and Jl^^ / Jl^^ = 0.12 for 7 = 20° 
leading to cj^ob — 3 tUrot, which is quite different from the experimental data. In contrast, 
as shown in Fig.0(b), the three moments of inertia calculated microscopically depend on 
cUrot even when the shape parameters are fixed, and the resultant cj^ob can either increase 
or decrease in general. In the present case of ^^^Lu in Fig.^ — J'^^^^ mainly determines 
the tUrot dependence. Its decrease is a consequence of that of J'x] the partial contribution 
to J7x from the proton ^13/2, ix/^rot, decreases as c^rot increases since this orbital is already 
fully aligned and therefore the aligned angular momentum is approximately constant. 
Thus, our result for cj^ob stays almost constant against u^ot, and even decreases slightly 
at higer frequencies approaching the experimentally observed one. This clearly shows that 
microscopic calculation of the three moments of inertia is crucial to understand the Wrot 
dependence of t^^wob in ^^^Lu. Let us compare this result with that for ^^"^Gd above. In ^^'^Gd, 
Ji^^^Jx ^ 1, J'i'''^^ ~ and IQi'^/Qi'^l < 1 at 7 < 20°. The last quantity measures the 
rotational if-mixing. This indicates that this solution is essentially similar to the gamma 
vibration in an axially symmetric nucleus as mentioned above. In contrast, the result that 
Jj^^'^/Jx = 0.90, Jt^^lJx = 0.19 and IQ^VQ^^I = 0.78 for ^^^Lu at hoorot = 0.3 MeV, for 
example, indicates that this solution is more like a wobbling motion of a triaxial body. The 
wobbling angle shown in Fig.^(a) is 19° - 13° for the calculated range. It is evident that the 
present small-amplitude approximation holds better at high spins. We confirmed that this 
wobbling solution disappeared as 7 decreased. Another feature distinct from the gamma 
vibration is that the present solution exists even at An = Ap = 0, whereas it is well known 
that the pairing field is indispensable for the existence of low-lying shape vibrations. This 
is related to such a tendency that the moments of inertia approach the rigid ones, J'x > Jy 
for 7 > 0, as the pairing gap decreases even without aligned QPs. 

A significant point of the data in Refs. [0, |^ is that the interband EM transition rates 
connecting the states J (TSD2) to / — 1 (TSDl) were precisely measured. In Fig.^, we com- 
pare our numerical results with the measured ones in a form similar to those in Refs. |0, § . 
Calculated values for / (TSD2) ^ / + 1 (TSDl) are also included in order to show the 



staggering behavior characteristic to this kind of transitions |T2[. Figure ^(a) presents the 
relative B{E2). The data indicate huge collectivity of the interband B{E2), such as 170 
Weisskopf unit. Although the present RPA solution is extremely collective, |c„=„ob| — 0.9 in 
the sum rule (Eq.(4.30) in Ref. [1^), in comparison to usual low-lying vibrations, the calcu- 
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lation accounts for 1/2 - 1/3 of the measured strength. Figure H(b) graphs B{Ml)/B{E2)ia. 
The smallness of B[M1) also reflects collectivity, that is, the coherence with respect to the 
E2 operator, indirectly. After confirming the insensitivity to g^^^\ we adopted 0.6 g^^^^"^^ 
conforming to Ref. |^ and calculated B{M1). The result is similar to that of the PRM. We 
confirmed that the sign of the E2/M1 mixing ratios was correct. 

To summarize, we have performed, for the first time, the RPA calculation in the ro- 
tating frame to the triaxial superdeformed odd-A nucleus ^^^Lu and discussed the physical 
conditions for the appearance of the wobbling solution in the RPA. We have confirmed 
that the proton ii3/2 alignment is indispensable for the appearance of the wobbling mode 
in this nucleus. The appearance of the wobbling mode requires J'^ > J^f^'' [y^ J^"^''^-, but 
the moments of inertia of the even-even core exhibit irrotational-like gamma dependence 
and therefore can not fulfill this condition for positive-gamma shapes. Consequently the 
alignment effect which increases Jx is necessary. Quenching of the pairing correlation also 
cooperates with the alignment effect for making the gamma dependence rigid-like. 

We acknowledge useful discussions with I. Hamamoto and H. Madokoro. This work was 
supported in part by the Grant-in-Aid for scientific research from the Japan Ministry of 
Education, Science and Culture (No. 13640281 and 12640281). 
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FIG. 1: (a) Moments of inertia and (b) wobbling frequency (left scale) and wobbling angle (right) 
in the five quasiparticle state in ^^'''Gd calculated as functions of 7 at ^rot = 0-3 MeV. The dip 
around 7 = 55° stems from a weak fragmentation of collectivity. Note that the present method of 
calculation does not apply to 7 ~ 0. 
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FIG. 2: (a) Wobbling frequency (left scale) and wobbling angle (right) and (b) moments of inertia 
in the TSD2 band in ^^^Lu as functions of htOxot- Here the latter were given by normalized to 
J'x{cil/2) = 99.2?i^/MeV. The proton BC crossing occurs at huj^ot ^ 0.55 MeV in the calculation. 
Experimental values were calculated from the energy levels in Refs. ||^, ^. 
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FIG. 3: Interband transition rates for / (TSD2) — > / it 1 (TSDl) transitions as functions of 2x 
spin /, (a) E2 and (b) Ml. They are divided by the in-band E2{I -^1 — 2) transition rates. 
Experimental values were taken from Ref. Q. Noting that the states / + 1 (TSDl) are slightly 
higher in energy than / (TSD2) at / > 51/2 and B{Tx;I ^ / + 1) ~ B{Tx;I + 1 ^ /) at high 
spins, we plotted those for I ^ I + 1 at the places with the abscissae / + 1 in order to show clearly 
their characteristic staggering behavior. 
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